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ABSTRACT 


For many years the cylindrical dipole and circular loop have been 
thought of as the basic radiating elements. Due to its relative 
Simplicity the study of the cylindrical dipole is quite complete, but 
owning to the degree of complexity, progress in the theory of the loop 
has been much slower. In recent years the isolated circular loop 
antenna of arbitrary size has been the subject of mathematical and 
experimental analysis [7], [6]. It was not until 1970 that the first 
complete and general mathematical solution for the mutual admittance 
of two arbitrary size coaxial circular loop antennas was published. 
The problem was formulated by Adler in terms of coupled integral equa- | 
tions for the current distributions on the loops [1]. The equations 
were solved using the Fourier series representation of the current 
distributions. The study presented here is the next logical step 
beyond the work of Adler. Using his current equations the current 
distributions are calculated at each point around the two loops. The 
input admittance, electric field and power gains are calculated. 
Radiation patterns and the loop current distributions are milonten for 
two cases: (1) both loops driven and (2) one loop acting as a short 


circuited parasite. 
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I. INTRODUCTION 


A. HISTORY OF WORK ON THE LOOP ANTENNA PROBLEM. 

In 1897 Pocklington excited a thin loop by a plane wave and 
obtained a solution for the induced current in the form of a Fourier 
Series expansion [9]. Hallen continued the work by producing an 
expression for the current distribution in the driven case [3]. Due 
to singularities arising in the Fourier series coefficients, Hallen's 
technique applied only to electrically small loops. Storer attempted 
to surmount the limitations of Hallen's work and derive an expression 
for the current in the case of a large loop [10]. He expressed the 
series in integral form and evaluated the integral in the sense of the 
Cauchy princival value. Experiments of Kennedy showed Storer's cal- 
culations produced an unacceptable deviation between calculated and 
measured susceptance [6]. Wu reviewed the work of Hallen and Storer 
and discovered that the shortcomings arose from over-simplifying the 
problem [11]. That is they assumed the electric field intensity on 
the surface of the conductor was uniform around the circumference with a 
filamentary current flowing at the center of the conductor. 

Wu then modified Storer's solution by changing the one dimensional 
integral equations to two dimensional integral equations in order to 
account for the electric field variation around the circumference of 
the conductor. Wu assumed a longitudinal current flow on the surface 
of the conductor and allowed no circumferential component. He showed 
this to be equivalent to a uniformly distributed current density on the 


cross section of the conductor, which allowed the current to be treated 





as filamentary and concentrated at the center of the conductor. After 
this simplification the equations were again one dimensional but no 
longer required uniform electric field around the circumference. King 
applied Wu's modification to calculate the admittance of a single loop 
[7]. This proved Storer's work to be valid but only for the case of 
small loops. 

‘ The first attempt at solving the multi-loop problem was by Iizuka 
for the special case of two identical loops [5]. In 1970 Adler advanced 
King's work by devising a completely general mathematical solution to 

the mutual admittance of large arbitrary size coaxial loops. Adler also 


measured the admittance in the two-loop case and found close agreement 


with theoretical values [1]. 


B. SCOPE AND LIMITATIONS OF THE STUDY 

This paper completes the study of large arbitrary size coaxial loop 
antenna arrays. Current scemaiens and radiation patterns are 
calculated and plotted for various array geometries. Although only the 
two-loop situation is presented herein, the theory applies to the 
n-loop situation as well. Since the @ (polarized) mode electric field 
component is of primary interest, this mode is examined in detail. The 
@ mode is calculated and one representative field pattern is plotted. 
The range of loop dimensions selected is for the resonant and near 
resonant case. a, b, and c are the conductor radius, loop radius and 
loop spacing respectively; the range of parameters are expressed as 

Wee s Kb = 2.0, OnS=< ke =< 720 


where: 


eel 
> 





ka 1S determined by the thickness parameter: 


2 = 2 In 2lb and is the wave length. 
a 


Having 7.5 s 2 < 15 represents a ratio of b/a from approximately 10 to 
300. The above parameter limitations are a result of the approximations 
used in the theoretical analysis of the Fourier coefficients of the 
current distributions. 

Although perfect conductors in free space are specified, the 
effect of lossy media can be included merely by considering the propaga- 
tion constant complex. The loops are driven by delta function voltage 
generators, symbolized by V 6(¢). The two coupled antennas are treated 
as a two-port network and conventional circuit analysis techniques are 


employed to produce a convenient matrix of equations. 





II, THEORY 


The conductor radius is assumed to be small relative to the loop 
radius, which allows the current density in each loop, 1, (¢) and I5($), 
to be approximated by an equivalent one-dimensional filamentary current 
distribution at the center of the loop. The electric field is derived 
by the scalar and vector potential in terms of the distribution of the 
current. Readers are referred to Adler's work for the mathematical 
details of the current equation formulation. Only an outline of the 
solution is presented here. 

As shown in Figure 1, the cylindrical coordinates for points on 
loops 1 and 2 are (by, $',0) and (b5,¢,c) respectively. The angle ¢' 
represents the coordinate of the current element, while the angle ¢ 
represents the electric field observation point coordinate. The 
distance from the source current element on one loop to the electric 


field observation point on the other loop is: 


R12 = Ro] = W (by -bo)* ete cf cts 4b, bosiné (26° : 


The distance from one point to another on loop 1 is: 
Roy = 4b- sin2 (oro) + 4a2 sin? y 
Vt ] 2 ll ‘ 


as shown in Figure 2. The same expression applies to loop 2 if the 
subscript 1 is changed to 2. yw is the angular distance around the cir- 
cumference of the conductor surface. In order to simply the notation 
in this development, the notation for harmonic time variation, eJut 


has been suppressed. 








Figure 1. Geometry of Two Unequal Loops 








Are at ya dws bj 


Figure 2. Geometry of Points on a Loop 
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e is the total electric field at a point due to loop currents 
1, (¢) and I5(>). The boundary conditions on e at the conductor 
Surface require the tangential component Ey to be zero at the surface; 
Since no circumferential current is assumed, the E. component is the 
only tangential component present. The delta function driving 


generator placed in the loop causes a singularity in the electric field 


at the gap such that: 


II 


6 in VRE (>) (1) 


| 


and hence: 


The boundary conditions on both loops are analogous. ec, is the ¢ 


> 
component of the total electric field due to the currents on both loops. 


These boundary conditions, when applied to the expressions for the total 
electric field intensity, specify the following coupled complex integral 


equations for the unknown current distributions: 


V, (4) = wre of [K,, (o-o')1, (o') + Kyo lo-o')Io(o') Ido: (2) 


V56(o) = =~ =f [K,,(o-¢ My (>') + Koolo-o' )Io(o') Ido! (3) 


In these equations 


—_ Pei. eee - 4! 4 
Kij(o-o') = [kbjb; cos(¢-9') + 7 a8 JW; 5(6-9') (4) 


2 





where Se is the kernel of the integral equations and 


[kb.b. cos(¢-9¢') + Lary 
j K 442 
p 
operates on Waa (o-9'). Wu's contribution to the solution resulted in 


integrating the term: 


ei to get: W..(9-9') =~ _— d 5 
where Wi; (o-o') is the Fourier coefficient of: 
aR 
a When i # j 
11 
SakKh. 
W..(o-¢') =F.) (6) 
J Rs j 


No is the characteristic impedence of free space. 

Equations (2) and (3) are now solved for I, (9) and 1,(¢) by the 
orthogonal expansion technique. The kernels of the integral equations 
and the currents are specified in terms of their respective Fourier 
Series as follows: | 

-jng' 


I.(¢) = qo le (7) 


age ij -jn($-9') (8) 
Ki jbee ae were 


n=-© 
By applying equations (7) and (8) to equations (2) and (3), inter- 
changing the order of integration and summation and performing the 


integration over ¢, equations (9) and (10) are formed: 





) i Als Ing 
eto, Ita, Ide = ae) (9) 


Np 


Meee sell <1. wee 
jx- } [ce 1 + Geet de V,6(¢) (10) 


After multiplying both equations by e INP and integrating over $ from 


-I_ to 1, the following simultaneous equations are obtained: 











1] 12 ] 
a a I V 
n n n _ ] (1 ] ) 
21 22 (2 a 7 ie 
=n on n 2 
Equations (11) are solved for ie and 1¢ and substituted back into 
equation (7) to get: 
12 
V, Gy 
22 
re Ying 
1,(¢) = a in. = ee (12) 
unio? if 12 
On OH 
aa! acc 
and 
: afl Vy | -ang (13) 
iia) = = TCOsmnG —ee— Te , 
Z JIng 7 
gle 
ra a2? 


The task now is to evaluate the ad's. The coefficients ald are put in 


terms of the Fourier coefficients <M of Wi 3(o-9') by letting: 


14 





aks a Inlo-d ) 


Ws s(o-9') = L : (14) 
where: 
ee oe Ve ee.) ade) 
Kn -_ — he o')e >. (15) 


The following difference relationship is arrived at after differentia- 


tion of equation (15): 


fee ig, KOGDS | 19 | | eB 
mn 1 aus es y 1 k “n - (16) 
a) = kb;b, «7, (17) 


For the evaluation of an refer to appendix A. 


The coordinate system is as indicated in Figure 3. The expressions 
for the far-zone electric field and gain are developed around Adler's 


coupled current equations (12) and (13) rewritten below: 





























2 12 
V a Vy a 
= ] 22 22 
I(¢) =-3 V5 o © Vo (18) 
Jing <i + 2) cos(n¢) 6 ee 
oo aoe “AO, n=] a ac? = oral 
1n() = LE] a2! vo c a2l Vp .(19) 
Jing | ———_——_—__—_- + 2) cos(ng) 55 TT 
all aoe -o14 2 ac! n=] al! af ~o,h2 On 


Let 1,(¢), the current in loop 1 conductor element b do' at the angle ¢' 


be as given in (18) above and I5($) be the same except for loop 2. 








(C) 





x 


Figure 3. Vector potential for one loop due to element of 
current I(9') b do' (a) at arbitrary point and (b) 
at point in plane of loop; (c) spherical coordinates. 
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The element of vector potential dA, at the point (x,y,z) due to 


1, ($') is composed of two components, dA,, and dAyy where: 








dA, = -dA, sin 9! (20) 
and 
dAyy = +dA, cos 9'. (21) 
So that: 
a7 JKR 
dA, , == HG ) z— Sin ¢ b1d¢ (22) 
and 
U aS 
= ' i ‘ 
where 
R= W2%+5% = (r, cos 9)¢ + s¢ (24) 
and 
s¢ = ce sin 9)¢ + be - 2br, sin 8 cos (p-o') . (25) 


When pecer ae ee Deron) COS (o-o') . 


The vector potential from loop | at points where bé<crg” is thus given 








by: 
bf .-Jkr 21 
= . a ec 2 (yj KBemS Ince Cos! }=¢." ). ae 
ALY ae Ju No aa sin o'dgo' (26) 
bu -jkr fii 
ae ociee =O 1, ,dkb, sin 6 cos($-9') ‘de! . (27 
er al ro F 11 (¢ Sie | CoS ¢da. = em) 


The same procedure is followed to get the vector potential radiated 


by loop 2. 
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Equation (18) is substituted into equations (26) and (27). The 


resulting integrals are solved as follows: 








12 
V5 any 
el 22 
1 (0) [Io of (28) 
Di, 


where I, (0) is the nth term in the series for loop 1 current at 9=0 


n 
(the feed point) and c =/“u/e where e is permittivity. 


DA = at! ane -al? ad! , (29) 
b,u  -jkr 
_ by ecdkrg 
i ee eee (30) 
B, = kb) sine . (31) 


Express sing’ and cos¢' in exponential form so that 


Q; 
Py a | 32 
ly j2 ( a b) ( ) 
n 
and 
Are = “Ty (K. + K,) (33) 
on a b 
n 
where 3 
2m j(n+1)¢' jB,cos(¢-9") 
af e e dq' (34) 
0 


1s j(n-1)9" jB,cos($-9') 
e a 


Ky = do' 5 (35) 


0 


These integrals are of the form of the Sommerfeld representation of the 
Bessel function: 


of" jBcos® jn¢ 
ae) 2m. e e do. (36) 





Letting ® = ¢-9' gives: 
--n ~jnd r2i-¢ jBcos(¢-9') jng' 
J, (B) -Je f e e 


do’ . 
2II -$ y 


This requires that 


n+] j(nt1)¢ (By) 
K. = 2 j é J+] : 


and 
Te Signin |) (By ) 
Kp = 2] J e Jn-1 


The components of the vector potential due to the component se 


Spherical coordinates are: 


Ay = (Ay cos ¢ + A} sin 9) cos 6 
0 an Yn 
n 
and 
AY = (-A] sing +A} cos 9). 
‘, Xn Yn 


Substituting (30), (34), and (35) into (40) and (41) gives: 





r JQ, : : 
7. sae [(K, -K,) cos¢ gS +Kp)sing] cose 
JQ] 
= = aoe ~Kyed?] COS¢ 
AD os In 6¢K +K,_) coso -j(K, -Ki,) sing | 
lp, la: ab) eae INNA O Sb > 


Q : 
= __f -J¢ Jo 
5 (Ke +Kpe ) : 


Substitute (38) and (39) into (43) and (45) to get: 


jn( ot as |e a paey (By) cose 


ee 
a? = -] Q, e n-1 


-2 Qy ein(or/2) n_ 3 (81? cose 
n B, 


and 


us 


(37) 





) in( +H 

“a a. a — ie, PP 7 Fe ees 48) 
| in($+I/2) 41 (B 

oe gy cteete) yO) (49) 


By summing all of the contributing vector potential components and 


uSing the recurrence relationships: 


2 2 
20h? = J (x) ft i (x) 


n+] ] 


and 
: (x) (x) 
Ae |) og ° 


the vector potential at a point due to loop | is: 














coe) foe) jn(o+ll B 
Ay = } Ay = -2m } a, e : Ban! )) cose (50) 
o -= 6, -0o ] 
and 
00 00 jn(o+m/2) B 
A, = } Al, = -jel U1 c a 1) (51) 
d =o n =O 
_ b, -jkr, 
Substitute Q, = =ly_ te —_~ into equations (50) and (51). 
4i¢z Yo 
V at? 
' = ; inlet) | J (kb, sine) 
‘ al -2Tl Qy : Vo On e kby Se cosé ° (52) 
D 
n 
Vy al? 
AY = -J20 Q, J [V2 of ]d! (kb, sine) . (53) 
> ac Dn 


Since a!' = a!! and J_,(*) = (-1)9 J 


=, n (x) | equations (52) and (53) can 


be rewritten as: 


20 





Vy elie 























n 
n 

AY = -21Q, Jn {V, a22 
1 1 2 “ni! J,(kb] sing) ar 

D, kb} sine 

jn(g+m/2) (-1)" -jn(otn/2)_ 
-[e - c ] (54) 
and 
12 
*, = -J21 Q)2 [Vp 082] Ji(kb, sine) + J [V5 af@/ 94 (kb, sine) 
or oa 1 

D. D, 


jn(o+m/2) -jn(o+m/2 
+( e ] 


e =1)" (55) 
The electromagnetic field from loop | is then given by: 
os ine Ayee eA) (56) 
] 9 Vy 


King shows that for a single loop of moderate size, |kb|<2.5, the 
Significant contributions are from the first three terms of the series 
[7]. 

The electric field patterns for the and 9 modes for loop | are 
readily obtained using the first three terms of the respective series. 
The identities: 


: 33%) z: = 5) 


and 





are used. 


21 




















n=0: V) 2 a}2 
a ree 22 
le 0; Mi, = j2m Q, [V2 26°] J, (kb, sine) (57) 
Do 
Vy a2 
ieee ee K : 
mig * Fl, * HO 0 2 a2] J1(kb,sine) (58) 
Do 


(Note that for the constant current distribution situation the series 


has only the n=0 term.) 



































n=1: Vy a4? 
AY = 41 Q) [Vo Vp af® J) (kbysing) cose sing | (59) 
S) —_—_———_—— ee 
] kb sine 
D, 
Vy a | 
,. = 41 Q Vp aie 22| Jy (kb, sine) COS¢ (60) 
D, 
r _-cr “ 
cilia 7° le, + 16, (61) 
12 
Vy 2) 
EY, =-J4m Qyu [Vp af@| Yi(kbisin®) cose sing (62) 
8 re or iS a 
] kb, sine 
D, 
12 
Vy Oy 
a me 254 V, af2! J! (kb, sine) (63) 
ly, = -J4 Qquw [V5 a} 1 Sin) cos¢ 
“4 
n=2: Vy aS* 
Ar = 34m Q, [Vp a5? 2Jo(kby sine) cose sin2s (64) 
89 — kb, sine | 
Do 


Ac 


























IZ 
V, a5 
»" jan Qe V. a5| Js (kby sine) cos2q 
Do 
ey =eES + 6 EL, 
2 9 2 
12 
Vo 
a = (4m Qqw) | V5 aS? 2Jo(kbySiN@) cose sin2¢ 
2 kby sine 
Do 
12 
Ef = (41 Qiu) |V5 aS¢| JS (kb, sine) cos2 
ly, 1 2 %9 2 1 > 
Ds 


The field from loop 2 is found in the same manner. 


(65) 


(66) 


(67) 


(68) 


To find the 


total field at a point for each mode, each term of the series for loop 2 


for both modes is multiplied by the phase factor, exp (jkc cose) to 


account for the separation of the loops. 


the » mode and @ mode are expressed as follows: 


and 


jkc cose jkc cose 
ce, - Ei, +65, e +E}, + ES e 
Pq = “bo °y ?] 
jkc cose 
+ Ey 7 E5 = 
°2 $2 
jkc cose jke cose 
ef = E£} Ep e + Ei, + E> e 
89 8 1 0] 
jkc cose 


r r 
+E +E e 
Ye, © “265 


Go 


The total electric fields for 





The gain of the array is calculated by dividing the square of the field 
by the power radiated. Where: 
. _ ae. cae 
Power Radiated = P44 = V4 YIN, + V5 YIN : (69) 


In order to find the input admittance, it is necessary to first derive 
expressions for self and mutual admittance. By treating the two coupled 
loops like the two-port network in Figure 4, equations (70) and (71) can 
be written directly: 


(@) © eS aes (70) 


1,(0) = Yoq Vy + Yoo Vp (71) 





Figure 4. Two-Port Network 


I, (0) and 1,(0) are the feedpoint currents. Because the network is 


reciprocal, Yy5 = Yoq° Solving for mutual admittance: 























12 ie 
: I, (0) _ ] “Ao ~ on? 
"2 ~ =a || ee ae 
V> Vj=0 JN 0 n=1 “n Voz 
Self-admittance is by definition: 
io) 
ian l (73) 
Vy V5=0 
I,(o) 
hee V Ve (74) 
Z 1 eames 
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The input admittance is by definition: 


I V 
Di eri 
s , ly 
Yay = eo HZ OY 
(ae SaieetS 
Pe Ws Vy 


(75) 


(76) 


In order to simplify the digital computer program used later in solving 


for the electric field and gain, it will be assumed that both driving 


voltages are purely real. This will allow YON to be replaced by Gyy 


where Gin equals the real part of Yr: 





Then: 
a 2 
Gin = Gy], + — 4&2 
1 Vy 
and 
4 
G =G,, + —G 
IN. 22 Vo ZA 
Finally, 
Gain, = |£), + esl’ 
Pad 
and: 
, r ry! 
Gain. = a eee 
’ Fi, * F2,| 
rad 


(77) 


(78) 


(79) 


(80) 


ED, and ED represent E2, and ED respectively multiplied by the phase 


+jkc cose 


factor (e ) which accounts for effects of loop separation. 


a, 





ITT. DISCUSSION CF NUMERICAL TECHNIQUES AND COMPUTER PROGRAM 


The ald terms were computed on an IBM System/360 Model 6/7 digital 
computer in FORTRAN IV language using Adler's basic program [1]. Beyond 
this point the programming was straight forward. Flow charts on 
Figures 5a, 5b, and 5c show the major steps in the program. The flow 
chart of Figures 5a and 5b is for computing I,(¢), I5($); Venice (up 
Y99s and Y12 where Meg and ee are the self admittance of loops | and 2 
respectively with mutual coupling effects ignored. Figures 5a and 5c 
show the steps in computing Gain, and Gaing. The computer results for 
current in the isolated loop were compared with King's [8], while the 
gain was compared with the work of Harrington and Mautz [4]. In both 
cases the agreement was exact. The incrementing of parameters b and c 
in the program allows for changing loop size and loop spacing. 


Incrementing 9 and ¢ allow for computation of gain patterns. 
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mee. -Ke , K 
fe) O O 












CALCULATE 
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n n 
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Figure 5a. Flow Chart of Computer Program 
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Figure 5b. Flow Chart of Computer Program 
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Figure 5b. Flow Chart of Computer Program (Continued) 
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Figure 5c. Flow Chart of Computer Program 
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Figure 5c. Flow Chart of Computer Program (Continued) 





IV. DISCUSSION OF RESULTS 


A. CURRENT DISTRIBUTION AND GAIN PATTERN CORRELATION 


Figures 6 and 7 show the patterns of Gain, in the x-z plane and the 


) 
X-y plane respectively. In all cases presented kb,=1 and loop 1 is in 
the x-y plane. Loop 2 is positioned coaxially above loop 1 at a 
distance c. Figures 8, 9, 10, 11 and 12 show the corresponding current 
distributions for the case of both loops driven in phase by one volt 
delta function generators. 

For .5 < kc < 1 and .5 < kbp < 1.5 the radiation patterns do not 
change significantly and only a slight change in phase of loop 2 current 
is noted. When kc > 1 radical changes in pattern shapes occur. Tnis 
Shows that as the separation becomes larger than about A/6, the effect 
of the array factor predominates. When kby = .5 the circumference of 
loop 2 is A/2. In this case the loop supports a 1/2 cycle standing 
current wave. 

Figures 13 and 14 are the Gain, patterns in the X-zZ and x-y planes 
respectively, for the case when loop 2 is parasitic. Figures 15, 16, 
17, and 18 are the corresponding plots of the current distributions. 
When the two loops are closely coupled, ke = .5 and resonant the 
currents are large and almost equal in magnitude but 180 degrees out of 
phase. Even though Gaing is relatively high in this configuration, the 
benefit of large currents is not realized because of the phase differ- 
ence between the two currents. As the separation of the loops 15S 
increased from approximately A/13 to approximately A/6, the currents 


drop by a factor of about five; however, the Gain, drops only slightly. 
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Sz 





Kc 5 A 72 6.3 


kb 


Zs 


‘Hat 


Figure 6. Gain ¢ patterns in x-z plane (= 0°,180°); feed point 
at g-0 , 6-90 3)¢-0 at Vigna. e-0 at top; 8-10, kbj=1, 
Vj=V0=1; kc=.5,1,3,6.3; kbo = .5,1,1.5,2. 
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Figure 7. Gain » patterns in X-Y plane, 9=90°; feed point at 
bottom (¢=0°); 2=10, alee gels kc — 0.153, 
6.3; Kis = 2 oainroace 





OvT 


O9T OT O21 001 


"Gg" =Cqy (9) 


AO} SUOLINGLUYSLP QUBAUND Z pue | sdoo, uo} aseud pue apnqLubey 


08 


09 0” 





02 


OT od 
O1Z- 


O8t- 
OST- 
OZT- 


06- 


¢ 
@ 


t=qy (q) tg°=%qy (2) § 


seeitep ((%1) eseuz 


881i 0desy{ 138 |< 1] 


Ost 


(q) 


O91 OFT O71 O01 


o] 


eal 


6 
0g 09 OF 


cmt 


G°=04 ‘|= 


02 


L 


0 


qy * Ls 


077- 
O1Z- 


O8t- 


OST- 


OZT- 


06- 


of 


09 


sevitep ((%)1] esvqg 


sssedeesiise |(%)1| 


C 


nz=tn *gL=0 


"8 aunBbl4 


077° 
O1Z-° 


ost- 
Ost¢ 
02T- 


Ob6- 





O91 O91 Ovl OZ7T OOT O08 09 OF O02 O 


C 


a00step ((M1) evrw 


sezedeny117s ((9)1| 


So 


~~ 2 





‘C'|s 


AOf SUOLZNGLUISLP JUBAUND Z puke | sdoo, uo} aseyud pue apnzrubew 


(2) 


O8T O9% OT OFT OOT 


4 
08 


Z 


09 0% 


0 


0 


qy (9) ° 


097° 


O1z- 


O8t- 


ost- 


OzI- 


06- 


09- 


ot- 


oc 


09 


2 marx , 
1=q4 (Gq) $5=%qy (e) f1=94 Stabqy <peeqely “gy =u 


(q) 


Ovz- 
Olt- 


OR T- 
OSt- 


OzI- 


09- 


eeortep (()1) serqu 


ot- 


Vyceat o¢ 
i ] 
r 


O8t O9T OT O71 OOT O08 OF OF OF O 


eezedmettiy= |(*)1] 
NI 
vr 
= 
“ 


seer3ep ((%)1) eevwm 


eesodestyty= |(9)1] 


Oet 


"6 oundL4 


(e) 


O9T OFT OZT OOT O08 


09 09 


0% 





0 


Oh dg 
Oz- 


ORt- 
Ost- 
ozt- 
06- 
09- 


ot- 


o¢ 


09 


- seeatep (()1) ovrqy 


sosedmeytty= |(%)1] 


36 





“S*L=Cay (2) 


1=¢qy (q) 


s 


$o°=Cqy (e) °€=0> ‘rolay epaCpaly ‘SL =8 


AO} SUOLINGLUZSLP JuUauuNd Z pue | sdoo| Ao} aseyd pue apnyiubey “o| aunbly 


(>) 
Orz- 
olz- 
oet- 
Ost- 
ott- 


06- 


¢ 


O81 O91 OT O21 OOT OF 09 OF OF O 


(q) 


Ovz- 
Olt- 
ogt- 
ost- 
oz1- 
06- 


09- 


vaar3ep ((%)1)} esryg 


ot - 


oC 


—— 09 


om 


Og! O91 O71 071 OO1 08 09 OF OF O 


corodeeytTy= |((%)1| 


paei8ep ((%)1]} sveeus 


corodeeytyys |( 1] 


os! 


091 


Ov 071 


001 


(e) 


ce] 


6 
0g 


09 OY 


07 


0 


072° 
Olt- 
OeI- 


OSt- 


OtI-. 


06- 





— Sereda stirs ((%)1| 


ee018ep (C91) sevw 


a7 





ost 


o? 
O9T OFT OTT OCOT O08 


09 OY 


6 
* 


9) 4 


L = qx (q) 


SUOLINGLUASLP JUaU4ND Z pue | Ssdoo, voy aseud pue apnzLubey 


0 


Ovt- 
OlZ- 


ORt- 


OStT- 


O@T- 


06- 


09- 


Ot- 


oC 


09 


aeei8ep ((%)1) sevq 


caredeesy17@ |(%)1] 


¢ 


S° 


= gy (e) ‘12a *1=%n='y *g=04 *G1=8 405 


(q) 


O8t O9T OT OZT COT O08 


09 OY 


02 


0 


OT 
Ole 


O8t- 


OST- 


OzT- 


oC 


09 


oS Cay ey 


eovitep ((%)1) vovqy 


serodeusy13@ ((*)1| 


Ost 


o9t 


“LL aunbly 


Ort O21 


oot 


(e) 


Qo 


¢ 
03 


ee 


09 OY 


oz 


0 


Ov 


Olz- 


Ost- 


OSt- 


Ovt- 


06- 


09- 


ot- 


oC 


09 


38 


eesatap (()1) evra 


seredevst {se [(*)1| 





ost 


O9t O71 O21 


ool 


° 


¢ 
08 


SUOLINGLUYSLP JZUBUUND Z puUeR 


oF 07 


02 


> [= 


¢ 


OWt- 
OIt- 


O8I- 
O¢I- 
O?I- 
06- 


09- 


Of- 


0 


oc 


09 


qy (q) 


eoer8ep ((%)1) eseqy 


evzedeeyiise |( 1] 


g's 


oo 


¢ 


7 


091 


(2) 


Ort O2T O01 


t) 


+ {=9) ‘|= 


Gal 


l 


gle | 


¢ 
0a 


09 07 


02 


0 


07%- 
O1t- 


O8I- 


O¢I- 


O?I- 


06- 


09° 


ot - 


o¢ 


09 


gy ‘L=¢A= 


eo0r8ep (()1} eerq 


cosedmestis= ((*1| 


¢ 


"6° t=°qy (9) 


Ly <GL=¥ aog 


O8t O9f Ont Oz1 


oot 


08 


| sdoo, wos} aseud pue apnz Lube, “2, aanbl4 


09 07 


é 





th) 
Olt 


OvI- 


O¢I- 


O@I- 


06- 


09- 


ot- 


oC 


09 


eeer3ep ((%) 1} oerw 


sosedemtttye |[(9)1) 


Sd 





Or 
Ase 


QO. 
=0° at top; =10, 


at =0°, 6=90°, o=0° at right, 6 
SD oa sos gl nonce 
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Figure 14. Gain patterns in X-Y plane, 9=90°; feed 
point at bottom (=0°; 2=10, Vy=1, Vo=0, 
kby=1; ke=.5,1,3.6.3; kb>=2551,1-5,2. 
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This observation is explained as follows: although the currents are much 
smaller when c = A/6, the phases of the two currents are no longer 
exactly opposing. This in turn allows the total vector potential of the 


¢ mode in both cases to be almost equal. 


B. CURRENT DISTRIBUTIONS 

Referring again to the current plots, the following observations 
are made. When both loops are resonant and driven, the current peak 
increases aS c increases from A/12 to A and decreases as c increases 
beyond }. When c is increased beyond 2A (not shown) the mutual coupling 
is negligible and the loop currents are very close to those of isolated 
loops. This phenomenon is due to the behavior of the input admittance, 


where: 


Yu, =Y,,+ 7 
IN; 11 12 Hy 1. 


By computing Yyq: Yoo: Yi0 and "INT 2 the following was noted: (a) with 
close coupling ‘10 has a large negative real part which approaches zero 
as c is increased; (b) the real part of Y11 is always positive; (a) 

Yry increases as c increases. Therefore, since 1(0) = Yyy¥, the current 
increases as Y;y increases. 

Of particular interest is the behavior of the current in the 
parasitic loop. In all cases both resonant and non-resonant the 
magnitude of the parasitic current, although reduced, paralleled the 
magnitude of the current on the driven resonant loop. In the resonant 
case when the spacing was close, the two currents were 180 degrees out 
of phase, but as the spacing approached i, the currents approached the 
in-phase condition. Also in the parasitic case the phases of both 


currents tend to be constant with an abrupt reversal at the nulls. 
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C. GAIN PATTERNS 

Figures 6 and 7, the Gaing patterns in the x-z plane and x-y plane 
respectively are again refereed to. Ine these figures 2-110, V)=Vo=! volt, 
Kb)=1, kc=.5, 1, 3 and 6.3 from left to right, while kbg=.5, 1, 1.5 and 
2 from top to bottom. In all patterns presented "+" marks the antenna 
center and "~~" marks unity gain. With kb, fixed, the patterns are 
seen to change shape as kc increases from .5 to 6.3. This is due in 
small measure to the change in mutual coupling and primarily to the 
array factor. The smallest Gain, me the Z ditmection occurs at Kb) = 
kbp = 1 and kc = 3 (c = A/2); the largest Gain, in the z direction 
occurs also with the resonant array but with kc = 6.3 (c = 2). In the 
former case the vector potentials were nearly 180 degrees out of phase 
and in the latter both vector potentials were in phase. In the extreme 
case where kb» = 2 and kc = 6.3, lobing occurs due to the combination of 
the array factor and tne diverse current distributions on the two loops. 

Gain, for the parasitic case 1s illustrated by Figures 13 and 14. 
The higher gains and more interesting patterns are attained in this 
configuration. As would be expected, less lobing now occurs and the 
parasitic element functions as a reflector. Of special interest are the 
kb. =) |-omekc = | and KDo = 2, kc = 6.3 patterns. In the first of 
these particular arrangements the loop array pattern resembles that of a 
log periodic antenna. Figure 2] shows the effect of changing kb» to 1.2 
and varying kc from .6 to 1.2. When kby = kc = 1.2 the antenna is uni- 
directional. In contrast to this configuration, note the resulting 
pattern when kby = 2 and kc = 6.3. The pattern is almost square with 


the gain across the top nearly uniform. Figure 14 shows a section of 


4/7 





these same patterns in the X-Y plane. These patterns closely resemble 
those of the familiar single-loop pattern in shape. 

The effect of changing the b/a ratio is illustrated in Figure 19. 
The change in 2 is from 7.5 to 15, which represents a b/a ratio change 
of about 10 to 300. Although the effect is measurable, it is concluded 
that the loop array is not sensitive to changes in 2. 

Figure 20 shows the consequence of ignoring the mutual coupling 
when calculating the gain pattern (the dotted curves are with coupling 
ignored). In Figure 20a (kc = .5 and kb, = kb, = 1) the coupling is 
very significant, but as seen in Figures 20b and 20c (kc = 6.3, kbo = 1 
and 1.5 respectively) the coupling has little effect. 

Sample patterns of Gaing are shown in Figure 22. Although 


Significant, this mode is certainly of secondary importance. 
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kb 5 2 
: 5 | 
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Figure 19. Gain ¢ patterns in X-Z plane ($=0°,180°) for 
V,=V5=kb =|, KC=6. 32a KD5—-oriielsoses O=7..5,15; 
¢=0° at right, 9=0° at op. 
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Figure 20. Gain » patterns in X-Z plane with (solid lines) 
and without (dotted lines) mutual coupling 
CTTCCtS a ie Me aKenaleeG-Gemdeerigit, 6-0 at 
top. 
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Figure 21. Gain-.¢ patterns in X-Z plane (¢=0°,180°) for 
= Opava te Vo=0, ko, =1, kb5=1. (ae ee e Ola. 
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Figure 22. Gain @ patterns for 9=10, kby=1 in (a) 9=45°, 
225° plane; (b) o=90°, 270°; (c) 9=90°, 270°. 
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V. SUGGESTIONS FOR FURTHER STUDY 


Since this study is of a theoretical nature an immediate require- 
ment 1s to confirm the results by experiment. First the current distri- 
butions and next the field intensity patterns for various arrangements 
Should be measured. Also of interest in the two loop array is the 
consequence of open circuiting the parasitic loop. 

The current and gain equations then should be solved for the three 
element case. Of interest would be the consequence of arranging three 


Or more elements to make a log periodic loop antenna. 
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APPENDIX A 


Evaluation of Fourier Coefficients Kd 
Case I: 7 = j 


Combining equations (5) and (15) results in: 


al P ] If Il 7 has (81) 
n 2 =p a ae 
(2m)¢ J, Ra 





By substituting A. =aZer sing. in equation (81) and interchanging the 


order of integration one arrives at the expression: 


be 2a. M (A.) 
a = a ee ee dA: (82) 
- 2 —. 
O 
where: 
; ] I jne -k ee 
M(A.) = — e e dé; n20 (83) 
n 1 2Tl 0 R. ; 
17 


Storer showed that for the case of A.<<b.: 





-> f [25,(x) + J Ipn(x) J dx (84) 


where a (x) is the mth order Lommel-Weber function and J (x) is the mth 


order cylindrical Bessel function of the first kind. 


oo 





Wu devised an approximation method to evaluate equation (84) based 
on the facts that (a) the integral representation of: 
2 1 
n 2\5 
(R= K2) 


b§ 


1 
where K(x) is the modified Bessel function of the second kind of zero 
order and the integral representation of biM CA.) are similar and (b) for 
A. small the difference between bs M,(A;) and the small argument approxi- 
mation for K, 1s independent of A;. The consequence of this is for 


moderate size n, 


1 no 2 Kb 
beeen) = = 2K, (A; =" k° } + In (1 - ~ : + 
1 





2kb; 
ZC en [25,,(x) + j Jo, (x) | dex (85) 
0 
n-] 1 
where C. =y-2 eo aa” In(4n) and y is Euler's constant. 


Because A. is small, the Ky term in equation (85) can be slightly 


modified with a negligible loss in accuracy to: 


2 nA. 2.2 
Lil tara || sole Pe _ Koby 
(tafe ) =e (Ft) bf) as 


j 


By employing equations (85) and (86) equation (82) can be integrated 


to get: : n 


17 ] Aca ie. (87) 
<n * 2b. 0 (F} ly | oe 


2kb. 
tf [Oq(x) + § Jpn (X)] | 





and: 
ii. 8b; ekb, 
Ko “opp> ln Fak af [Q9(x) = j Jan(x)] dx? - (88) 


Case II: ifj 





By combining equation (15) with: 


| 1 
= ‘ et =o" ae j 2 2q5 
Ra [4b jb. sin (oe) + (b; b;) +5 Cal 


and 
ie. - t = eal) 
ij (o-9') Raj 


equation (89) follows: 


I ] 

e e 2 2 do 2 2 ron 

~ jng.-jk [4b.b. sin peepee. (Cm2 
ae) i Se) eee i 9) 


n 20 My ean ene 251 

si [4bjb; sin $+ (bj-b,) + c], 

It is necessary to solve equation (89) by numerical methods and the 
reader is referred to Adler's original work for the details of this 


process [1]. 
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